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Abstract 

Using a simple example of biological macromolecules which are partitioned be- 
tween bulk solution and membrane, we investigate T.L. Hill's phenomenological nan- 
othermodynamics for small systems. By introducing a systems size dependent equilib- 
rium constant for the bulk-membrane partition, we obtain Hill's results on differential 
and integral chemical potentials fi and fi from computations based on standard Gibb- 
sian equilibrium statistical mechanics. It is shown that their difference can be under- 
stood from an equilibrium re-partitioning between bulk and membrane fractions upon 
a change in system's size; it is closely related to the system's fluctuations and inhomo- 
geneity. These results provide a better understanding of the nanothermodynamics and 
clarify its logical relation with the theory of statistical mechanics. 

1 Introduction 

With great advances in material preparations at nanoscale in recent years, there is a grow- 
ing interest in the thermodynamics of small systems 03 13- Nanothermodynamics, a phe- 
nomenological theory developed by T.L. Hill in the early 1960s for small equilibrium sys- 
tems, has become one of the major formalisms for quantitative treatment of equilibrium 
nanoscale materials jUSS. There are many experimental and computational work that have 
verified the theory, see for examples [H|6]|. 

The greatest strength of thermodynamics, being absolute and abstract, is often its weak- 
ness: It usually provides no molecular mechanism and insights into chemical and bio- 
chemical processes. The same dichotomy applies to Hill's nanothermodynamics: The 
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key concept in the theory is a difference between differential and an integral froms of 
many non-extensive thermodynamic quantities due to smallness of a system. For exam- 
ple, since Gibbs free energy G(N,p, T) is not linearly proportional to the systems size N, 
fi = dG/dN and fi = G/N are different. The difference is intuitively understood: It is due 
to contributions such as surface effect. 

In biophysical chemistry, biological macromolecules are often treated as having dis- 
crete conformational states, and interaction with ligands as stoichiometric binding. While 
this approach is only an approximation, it often provides much more powerful intuitive 
understanding of the physical processes that underly molecular thermodynamics. A case in 
point is J. A. Schellman's theory of three-component system with mixed solvents 0. While 
a thorough treatment of this problem involves perferential interaction coefficients, discrete 
stoichiometric binding model has provided greater intuition (8ll3 

The fundamental premises of a small system is that its thermodynamic properties are not 
strictly proportional to the system's size [4|. A system's size, however, can be represented 
by several different quantities: The number of molecules N in a system; the volume V of 
a system, or the total energy E of a system. In chemistry and biochemistry, concentration 
x = N/V is a widely used intensive quantity. We shall refer to system's size in term 
of its volume in the rest of the paper. In terms of V, the Helmholtz free energy F of a 
three-dimensional small system could have the form 

F = Vf(x, T) + a(x, T)V 2/3 + b(x, T) In V + ■ • • , (1) 

where the term in V 2 ^ 3, is associated with a surface energy. 

In the present work, we first consider the following simple problem: iV non-interactive 
macromolecule M are in a system which consists of bulk solution and membrane. (One 
could think of a vesicle with lipid membrane.) Each M molecule has three discrete con- 
formational states: 0, 1 in the bulk and 2 on the membrane. We assume the equilibrium 
constant between states and 1 is K% and between states and 2 is K 2 . K\ = pi/po and 
K-z = P2/P0 where pi is the probability of an M being in state i in equilibrium. K\ is 
independent of amount of membrane in the system; but K 2 is critically dependent on the 
ratio of membrane surface S to bulk volume V, thus the system's size. Assuming simple 
geometry with surface to volume ratio S/V = V' 1 ^ 3 and homogeneity in both bulk and 

'A difference between differential and integral forms of macromolecular interaction also appears in the 
theory of binding: Thermodynamic binding can go to zero while molecular interaction is the strongest |9|. 
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membrane phases, it is reasonable to assume that K 2 (V) oc V 1 ^ 3 . Thus the equilibrium 
partition function for a single M molecule, with the state as the refererence state: 

Z x = 1 + K x + K 2 (V) ->• 1 + K x (2) 

in the limit of V — > oo. The simple example can be easily generalized to more general 
macromolecular partition problem, and a key formula relating n and jl first appeared in flU, 
will be derived. 

A few remarks on Eq. © are in order: (i) Z\ looks like a binding polynomial, but it is 
actually independent of the concentration of the M. This is because we have assumed that 
theiV macromolecules are non-interactive. Hence, Z x is the "partition function" of a single 
M molecule. The partition function for the whole system is Zn = Zf. 

(ii) K 2 (V) = p 2 /po as a function of V is itself a complex problem in general. For our 
simple model, we shall completely neglect surface free energy of the membrane. Further- 
more, if one treats macromolecule-membrane association as simple Langmuir absorption, 
with N(p + pi)/V and Np 2 /S being the concentrations of the M in bulk and on mem- 
brane. Then one obtains p 2 /(po + p\) oc S/V. Since (p + p x ) = (1 + K x )p , one has 
P2/P0 oc S/V ~ V~ x l' i . For the present analysis, the key is that K 2 (V) is a function of V. 
The result in Eq. © given below does not depend on a specific functional form of K 2 (V). 

The novel idea of the present work is to introduce the systems size dependent equilib- 
rium constant(s), (i.e., more generally potential of mean force or conditional free energy). 
It is shown that if one starts with a mechanistic based statistical approach, following Gibbs, 
to small systems, then the small systems nano thermodynamics naturally emerges. 

One of the salient features of Hill's nano thermodynamics is its ensemble dependence: 
When a systems boundary becomes significant, how a statistical thermodynamic system is 
maintained at the boundary matters. In the past, dependence upon details of an ensemble 
at its boundary leads to the phenomenon of entropy -enthalpy compensation ifTOlfTTl . which 
can be understood from relaxing certain internal constraints. In Sec. I3.1l we show that this 
perspective is also fruitful in understanding the difference between differential and integral 
chemical potentials fx and jl. 
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2 Bulk-Membrane Partition: A Simple Example 



We now consider a simple system in which there are N non-interactive macromolecle M. 
Each M has three distinct states: states and 1 are in bulk solution and state 2 is associated 
with membrane. It is clear that the ratio between total membrane to total bulk volume 
is a function of the systems size V: Thus by simple geometric consideration K 2 (V) is 
volume dependent. The partition function for the small system with constant pressure p 
and temperature T is: 

Z{N,p,T) = (l + K 1 + K 2 {V)f, (3) 
in which V = V(N, p, T) is a function of JV@ Therefore, 

G(N,p,T) = -k B TlnZ{N,p,T), (4) 
fi(N,p,T) = Gi ~ N ^ T) = -k B T In (1 + K 1 + K 2 ) , (5) 

(N dG(N,p,T) . k B TK 2 (») p<r 

^ N ^ T) = ON 1 + K 1 + K 2 ' (6) 

where (i and p, are integral and differential chemcial potentials introduced by Hill H|. 

In the last term of Eq. ©, (d\n.K 2 /d\nN) is an intensive quantity, K 2 oc V^ 1 ^ 3 and 
1 + K\ + K 2 — > \ + K\. Hence it is on the order of K^ 1 / 3 . Thus, p and p, are the same for 
a macroscopic system when G(N,p, T) asymptotically becomes an "extensive" quantity: 
G oc N. 

Eq. © can be easily generalized into 

(AT T\ "(AT ^ / ( dp°(N^T) \ \ f dp(N,p,T) \ 

where 

p°{N,p.T) = -k B T\nK t (N,p,T) (8) 

and the average (■ ■ • ) is performed with respect to i = 0,1,2. Eq. © clearly shows why 
and how the small systems thermodynamics depends on the nature of ensemble; and why 
the difference p and p will be ensemble specific. Eq. © agrees with what obtained by Hill 
(HI, Part I, p. 30, Eq. 2-12): 



2 Mathematically more rigorous, one needs to solve the problem with self-consistency among Eqs. (01, (0]) 

md:V(N,p,T) = (dG/dp) NT . 
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3 Thermodynamic Hierarchy and Nanothermodynamics 

3.1 Perturbation and relaxation in thermodynamic systems 

The chemical potential fi is the increment in free energy of a system when its number of 
molecules increased by 1 lfT2l . Therefore, it can be considered as a perturbation to the 
thermodynamic system. The idea advanced in IfTOl , originally proposed by A. Ben-Nairn 
|fl"3l , is to view the response of a thermodynamic system to a perturbation as two parts: 

(i) Constraining the distribution among all the populations un-changed, what is the 
increase in free energy of the system due to the perturbation? 

(ii) Inevitably, however, the perturbation will cause a shift (i.e., re-distribution, re- 
partitioning) in the equilibrium among the populations. There is an additional free energy 
change associated with this re-organization process. 

This perturbation-relaxation view of thermodynamics of complex systems yielded a 

possible explanation for the so-called "entropy-enthalpy compensation" phenomena |[T3l 

[TUMUl: Note that the chemical potentials for all the subcomponents, fij, of an equilibrium 

system have to be the same. Thus, the free energy change associated with (ii) is essentially 

zero if the perturbation is infinitesimal. In contrast, the corresponding entropy and enthalpy 

of the subcomponents, s~ and hj (hj — Tsj = uA can be very different. Thus, a large 

□ 

contributions to entropy and enthalpy from the process in (ii) are expected^ 

Applying the same idea to thermodynamics of our simple example, a perturbation leads 
to a shift in equilibrium between the "bulk" and "membrane" fractions. More specifically, 
introducing an additional M molecule causes a change in the system's volume V under 
isobaric (i.e., constant p) condition. This change in V is global: It causes the change of K 2 
for each and every molecule M in the system. 

3.2 Spatial inhomogeneity 

"Size" is a geometric concept. The unique feature of nanothermodynamics of small systems 
is spatial inhomogeneity that partitions the particles in a system into different populations 
whose proportions scale differently with systems size. In Hill's theory, this partition is 



3 Both steps in this gedankenexperiment are carried out under an isothermal condition. They can not be 
realized in a laboratory. They are different from the adiabatic (isoentropic) and isothermal processes in the 
derivation of the fundamental equation of thermodynamics. 
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implicit: cluster growth naturally leads "surface" and "interior"; in our simple model this 
partition is explicitly assumed: states and 1 versus 2. 

With this geometry in mind, there will be two answers to the following question: What 
is the change in system's free energy if one introduces one additional particle into the 
system? Accordingly, answer (i) assumes that the proprotion is unchanged; and answer 
(ii) assumes there is a shift in equilibrium distributions between different populations. 

The answer (i) is precisely the macroscopic result, fi. It is intimately related to the 
free-energy perturbation theory |[T4l[T5l : 

AG°(A ^B) = -k B T In ( exp ) A> (10) 

in which the average is carried out in terms of the equilibrium distribution of the unper- 
turbed system. Note that Eq. (flOl) assumes that the phase space, over which the ensemble 
average (• • - )a is carried out, is the same before and after the perturbation. 

The answer (ii) however, contains exactly a contribution from small systems effect: The 
addition of a particle changes the size of the system, thus the phase space, on the order of 
1/N. This change causes a shift in the equilibrium between "surface population" to "bulk 
poplulation", which affects all iV particles in the system. This is given in Eq. ©. 

Realizing the difference between these two answers led Hill to introduce dN t and dj\f. 
The former corresponds to a change in the total number of particles, thus gives \i (answer 
ii); the latter changes the number of particles but keeps the equilibrium population propor- 
tion by "adding an independent copy of the small system", thus gives fi (answer i). 

3.3 Rapidly stirred biochemical reaction systems 

For a small biochemical reaction system which is rapidly stirred, the theory of Delbriick- 
Gillespie process [fT6l [T71 . with its probability distribution given by the chemical master 
equation and its trajectories follow Gillespie algorithm, predicts that its stationary proba- 
bility distribution, in the limit of large size, has the generic form 

f(x)<xexp(-V<p(x)), (11) 

where x is the concentration of biochemical species, V is system's volume, and <p(x) is a 
function of only intensive quantities. Applying Laplace's method for integrals [HHHl, this 
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result leads to the conclusion that for a rapidly stirred system without geometric inhomo- 
geneity, the partition function 

« Ve~ v ^ J dxe-^"^ x - x ^ 

where is the macroscopic concentration. Therefore, 

- k B T In Q(V,fi,T) = k B TV(f)(x*) + 0(\nV), (13) 

where 0(- ■ ■ ) is the mathematical symbol called Bachmann-Landau notation. It stands for 
"on the order of • • • ". Eq. CCl should be compared with Eq. ©: -k B T In Z(N, p, T) = 
— k B TVx* ln(l + Ki) + (V^ 2 / 3 ) where x* = N/V. For rapidly stirred nanoscale systems, 
the correction to the extensive term is on the order of logarithm of an extensive quantity: 
In V or In N. When there is a spatial partitioning, e.g., compartimentalization, however, 
terms on the order of fractional power of the extensive quantity arise. For system with 
bulk-membrane partition, this term is on the order of V 2 ^ 3 . 

A biological cell is a small thermodynamic system. Applying the present result to a 
collection of cells, one realizes that there are two macroscopic limits: (a) Classical bio- 
chemistry studies cell extracts with membrane removed; (6) permeable cells that preserve 
the ratio between membrane surface area to bulk volume. They correspond precisely the 
two ensembles of Hill's. 



3.4 Kinetics of re-partitioning 

All the above discussion has been exclusively on equilibrium thermodynamics. We suspect 
there is in fact a kinetic aspect of the thermodynamics of small systems. Nonequilibrium 
thermodynamics for heterogeneous systems and for nanoscale systems have become matu- 
ral subjects in recent years [|20ll2Tll22i . If one takes a kinetic perspective of the response of 
a small system upon perturbation IfTOfl , one sees that the conformation re-partitioning be- 
tween the bulk and surface could indeed be observed kinetically. This connection between 
thermodynamics and kinetics deserves further investigation ||23l 1241 |25i 
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